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SUMMARY 


The  reciprocal  equations  relating  the  lift  and  moment  acting  on  an 
airfoil  due  to  harmonic  oscillations  in  compressible  flow  to  the  indicial 
lift  and  moment  functions  for  a  pitching  airfoil  are  determined.  Two 
indicial  functions  are  required  to  describe  completely  the  lift  or 
moment  on  the  airfoil;  one  function  is  used  for  describing  the  effects 
due  to  the  angular  position  of  the  airfoil  whereas  the  second  function 
is  used  for  describing  the  effects  due  to  pitching  velocity. 

Calculations  for  the  indicial  lift  and  moment  functions  are  made 
for  a  pitching  airfoil  at  a  Mach  number  of  0.7  for  which  sufficient 
oscillatory  coefficients  were  available.  The  growth  of  lift  to  the 
steady-state  value  is  less  rapid  for  compressible  flow  than  for  incom¬ 
pressible  flow  for  both  indicial  functions  of  the  pitching  airfoil. 
Although  the  circulatory  lift  due  to  pitching  velocity  on  an  airfoil 
rotating  about  its  three-quarter-chord  point  is  nonexistent  at  a  Mach 
number  of  0,  for  a  Mach  number  of  0.7  a  time-dependent  lift  exists 
because  of  pitching  velocity  and  approaches  zero  in  approximately 
2  chords. 


INTRODUCTION 


In  the  study  of  transient  flows,  two  types  of  airfoil  motions  have 
special  significance  -  a  harmonically  oscillating  airfoil  and  an  air¬ 
foil  experiencing  a  sudden  change  in  angle  of  attack.  The  reciprocal 
relations  between  the  indicial  lift  and  the  lift  associated  with  a 
harmonically  oscillating  airfoil  for  both  sinking  and  pitching  motion 
are  given  in  reference  1  for  incompressible  flow.  In  reference  2  the 
reciprocal  relations  were  extended  to  apply  to  an  airfoil  suddenly 
acquiring  a  vertical  velocity  in  a  subsonic  compressible  flow,  and  the 
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indicial  lift  and  moment  functions  were  evaluated  for  a  Mach  number 

of  0.7*  The  present  paper  is  an  extension  of  reference  2  in  that  the.  1  ^ 

indicial  lift  and  moment  functions  for  a  pitching  airfoil  are  determined. 

The  reciprocal  equations  for  the  case  of  a  pitching  airfoil  are  indi¬ 
cated,  and  the  indicial  lift  and  moment  functions  are  computed  for  a 
Mach  number  of  0.7. 

In  a  recent  paper  by  Lomax,  Heaslet,  and  Sluder  (reference "3),j  a 
different,  method  for  determining  the  indicial  lift  and  moment  functions  ; 

is  given.  Although  parts  of  the  indicial  functions  can  be  determined 
readily,  the  solution  for  the  complete  indicial  functions  is  lengthy 
and  tedious,  and,  consequently,  numerical  results  are  given  only  for  a 
Mach  number  of  0.8  in  reference  3*  Part  of  the  solution  presented  herein 
for  a  Mach  number  of  0-7*  however,  was  determined  by  the  method  in  refer¬ 
ence's  and  was  compared  with  the  solution  determined  in  this  paper. 


SYMBOLS 


s 

0) 

V 

c 

x,z 

k 

L(s) 

Lq(s) 

M(s) 

Mq(s) 

P 


distance  traveled,  half-chords 

angular  frequency 

forward  velocity  of  airfoil 

•  V 

chord 

Cartesian  coordinates 

■  ,  1 

reduced- frequency  parameter  (ox:/ 2V) 

lift  per  unit  length  of  span,  positive  in  downward 

direction  - 

lift  per  unit  length  of  span  due  to  pitching  velocity, 
positive  in  downward  direction  » 

moment  per  unit  length  of  span  about  quarter -chord 
point,  positive  when  moment  tends  to  depress 
trailing  edge 

moment  per  unit  length  of  span  due  to  pitching  velocity 
about  the  quarter-chord  point,  positive  when  moment 
tends  to  depress  trailing  edge  1 

density  P 


/ 
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3 


% 

it 


v 


* 


<i 


normal  component  of  perturbation  velocity 


h 

e 

q 


ki(s) 


mi(s)' 


k]_  (s) 

q 


miq(s) 

Fc(k)  +  iGc(k) 


amplitude  of  vertical  displacement  of  airfoil,  half- 
chords 

angle  of  pitch,  positive  when  trailing  edge  is  lower 
than  leading  edge 

1 

rate  of  change  of  pitch  of  airfoil  with  respect  to 
distance  traveled  in  half-chords,  positive  when 
trailing  edge  is  falling,  with  reference  to  leading 
edge  (6c/2V) 

indicial  lift  function  for  an  airfoil  experiencing  a 
sudden  change  in  vertical  velocity 


indicial  moment  function  for  an  airfoil  experiencing 
a  sudden  change  in  vertical  velocity  where  moment 
is  taken  about  quarter-chord  position 

indicial  lift  function  for  an  airfoil  experiencing  a 
sudden  change  in  pitching  velocity  about  its 
leading  edge 

indicial  moment  function  for  an  airfoil  experiencing 
a  sudden  change  in  pitching  velocity  about  its 
leading  edge 

coefficient  of  complex  compressible  oscillatory  lift 
derivative 


M(k)  +  iN(k)  coefficient  of  complex  compressible  oscillatory 

moment  derivative  where  moment  is  taken  about 
quarter-chord  point 


Fcq(k)  +  iGcQ(k)  coefficient  of  complex  compressible  oscillatory  lift 

derivative  due  to  pitching  velocity  only 

Mq(k)  +  Nq(k)  coefficient  of  complex  compressible  oscillatory  moment 

derivative  due  to  pitching  velocity  only 


f(k)  =  Fc(k)  -  Fc(oo) 


m(k)  =  M(k)  -  M(°°) 

fq(k)  =  Vk)  -  Fcq(°°) 
mq(k)  =  Mq(k)  -  Mq(°°) 
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k 


in-phase  and  out-of-phase  lift  coefficients, 

respectively,  associated  with  translation  of  airfoil 

Mi,M2  ,  in-phase  and  out-of-phase  moment  coefficients, 

respectively,  about  quarter-chord  point  associated 
with  translation  of. airfoil 

ift-phase  and  out-of -phase  lift  coefficients, 

respectively,  associated  with  pitching  motion  of 
airfoil  about  its  leading  edge 

in-phase  and  out-of -phase  moment  coefficients, 

respectively,  about  quarter-chord  point  associated 
with  pitching  motion  of  airfoil  about  its  leading 
edge 

METHOD  OF  ANALYSIS 

Consider  the  perturbation. velocities  on  a  two-dimensional  wing 
(see  the  following  sketch)  for  the  case  of  combined  vertical  motion  and 
pitching  moti6n.  For  convenience  the  reference  axis  is  taken  as  the 
leading  edge  and  z  and  G  are  used  to  describe,  respectively,  the 
vertical  and  angular  positions  of  the  airfoil. 


z 


The  perturbation  velocity  (w  measured  positive  upward,  z  positive 
downward)  may  be  expressed  as 

w(x,t)  =  -(VG  +  z  +  x0) 

If  the  airfoil  is  considered  to  have  vertical  motion  only,  then 
the  perturbation  velocity  will  be  uniform  across  the  chord  and  will  be 
of  intensity  z.  If  the  case  of  an  airfoil  pitching  about  its  leading 
edge  is  considered,  however,  the  perturbation  velocity  will  be  composed 
of  two  parts:  a  uniform  part  of  intensity  VG  and  a  linearly  varying 
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part  of  intensity  x@.  When  the  lift  on  an  airfoil  having  arbitrary 
motions  is  determined,  therefore,  these  observations  mean  that  only  two 
types  of  perturbation  velocities  have  to  be  considered:  a  uniform  dis¬ 
tribution  which  applies  either  to  vertical  motion  or  to  angular  position 
alone  and  a  linearly  varying  distribution  which  applies  to  angular 
velocity  of  the  airfoil. 

When  unit  step  motions  of  the  airfoil  are  considered,  it  should  be 
apparent  that  one  indicial  function  is  sufficient  to  define  the  lift 
on  an  airfoil  suddenly  acquiring  a  vertical  velocity,  since  only  the 
uniform  perturbation  velocity  is  involved.  For  the  case  of  a  pitching 
airfoil,  however-,  two  indicial  functions  are  required:  one  corresponding 
to  the  uniform  perturbation  velocity  associated  with  angular  position 
and  the  other  corresponding  to  the  linear  variation  which  is  associated 
with  angular  velocity. 

The  lift  on  an  airfoil  suddenly  acquiring  a  vertical  velocity  is 
generally  given  in  terms  of  an  indicial  lift  function  kq(s)  by  means 
of  the  equation 


_  L(s)  =  -*PcV2  ||  kx(s)  (la) 

where 

dh  c  »  _  z 
ds  ~  2V  n  “  V 

This  equation  may  he  made  to  apply  to  an  airfoil  suddenly  acquiring  an 
angular  position^  because  of  the  similarity  that  exists  between  the  per- 
turbation  velocities  for  the  case  of  vertical  motion  and  for  the  case  of 
angular  position  alone.  The  lifts  for  these  two  cases  therefore  will 
be  equal  if  the  intensities  of  perturbation  velocity  are  equal,  that  is, 
if _  z  =  VS..  With  this  condition  and  equation  (la),  the  lift  on  an  air¬ 
foil  following  a  sudden  change  in  angular  position  is 


L(s)  =  -JtpcV^k^s) 


(lb) 


For  the  pitching  case  then,  the  only  indicial  lift  function  that  remains 
to  be  determined  is  the  function  associated  with  pitching  velocity  of 
thd  airfoil.  The  equation  for  lift  following  a  sudden ■ change  in  angular 
velocity  0,  where  this  angular  velocity  is,  for  convenience,  taken 
about  the  leading  edge,  may  be  written  similar  to  equation  (la)  or 


Lq(s) 


-2tfpcV2qki  (s) 
q. 


*■  f' 
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is  the  indicial  lift  function  associated 
an  angular  velocity  about  the  leading  edge. 

The  indicial  lift  and  moment  functions  due  to  pitching  velocity 
can  be  determined  with  the  aid  of  the  appropriate  reciprocal  equations 
and  flutter  coefficients  in  a  manner  similar  to  that  shown  for  the 
sinking  airfoil  in  reference  2.  The  available  values  for  the  lift  and 
moment  on  the  oscillating  airfoil  for  pitching  motion  at  subsonic  Mach 
numbers,  however,  have  been  obtained  numerically  only  for  the  total  lift 
and  moment)  that  is,  the  components  due  to  angular  position  and  pitching 
velocity  have  not  been  separated.  Fortunately,  the  component  due  to 
angular  position  may  be  subtracted  out  of  the  total  oscillatory  coeffi¬ 
cients  because  of  the  similarity  that  exists  between  the  lift  on  a 
sinking  airfoil  and  the  component  of  lift  associated  with  angular  posi¬ 
tion  of  a  pitching  airfoil.  The  available  data  for  an  airfoil  oscil¬ 
lating  harmonically  in  vertical  motion  therefore  are  used  in  the  reduc¬ 
tion  of  the  data  for  an  airfoil  oscillating  in  pitch  to  give  the  lift 
and  moment  associated  with  oscillatory  pitching  velocity  alone.  These 
data  can  then  be  used  to  determine  the  indicial  functions  associated 
with  pitching  velocity. 


where  q  =  ^  9  and  ki^(s) 
with  a  sudden  acquisition  of 


Reciprocal  Equations 


The  reciprocal  relations  given  by  equations  (8a)  and  (8b)  in  refer¬ 
ence  2  between  the  compressible  lift  on  an  airfoil  due  to  harmonic 
oscillations  and  the  indicial  lift  on  an  airfoil  experiencing  a  change 
in  angle  of  attack  due  to  a  sudden  acquisition  of  vertical  velocity  are 
as  follows: 


o  r°°  Fc(k)sin  ks  .  . 

ki(s )-|j  (®  >  0) 


(2a) 


and 


o  n00  G„(k)cos  ks 

kl(s)  =  Fc(0)  +  \  /  - - 7T—  dk  (s>0) 


(2b) 


where  Fc(k) •  and  Gc(k)  are  related,  respectively,  to  the  in-phase  and 
out-of -phase  lift  components  on  an  oscillating  airfoil  and  are  defined 
by  the  following  equation: 


L(s)  =  -itpcV2eiks 


(-ikh) 


(k)  +  iGc 
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V 

1 


For  the  case"  of  a  pitching  airfoil  the  component  of  lift  associated 
with  angular  position  and  its  corresponding  indicial  lift  function  are 
given  by  equations  (2)  and  (3),  where  the  term 

-ikh  =  0  (k) 

a  relationship  which  follows  from  the  fact  that  the  perturbation  veloci¬ 
ties  for  vertical  motion  are  equal  to  the  velocities  due  to  angular 
position  alone  when  z  =  V0. 


In  the  interpretation  of  the  expression  for  lift  given  by  equa¬ 
tion  (3),  it  should  be  noted  that  the  expression  for  the  motion  of  the 
airfoil  is  in  complex  form.  If,  for  example,  the  actual  angular  posi¬ 
tion  of  the  airfoil  is  denoted  by 

I.p.eeiks 


then  the  lift  is 


I.P.-npcV20eiks 


Fc(k)  +  iGc(k) 


3 


or 


TTpcV2©  Jl 


FcOO  +  GcdO  sin 


1/2 


ks  +  tan 


_!  Gc(k) 

FXkJ 


Thus,  the  lift  and  the  angular  displacement  have  the  same  frequency  and 
both  their  magnitude  and  phase  are  functions  of  k. 

The  reciprocal  equations  for  the  component  of  lift  due  to  pitching 
velocity  may  be  expressed  in  terms  of  a  harmonically  oscillating  airfoil 
and  the  indicial  lift  function  due  to  a  sudden  change  in  pitching 
velocity  kpg(s)  as  follows: 


and 


kl  (s)  =  w 
q  n 


00  F  (k)  sin  ks 


dk 


0 


(s  >0) 


(5a) 


kl  (s)  =  Fc  (°) 

q  q 


Gc^(k)  cos  ks 


(s  >0) 


(5b) 


0 


k 


dk 
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where  Fc  (k)  and  Gc  (k)  are  related,  respectively,  to  the  iri-phase 

and  out-of -phase  lift  components  on  an  oscillating  airfoil  due  to 
pitching  velocity  about  the  leading  edge  and  are  defined  by  the  equation 

L  (s)  =  -7tpcV2eiks(2ik0)  Fc  (k)  +  iGc  (k)l  (6) 

”  .  L  uq  q  J 

In  equations  (5a)  and  (5b)  the  functions  FCq(k)  and  GCq(k)/k  must 
be  continuous  and  finite  in  the  interval  from  k  =  0  to  k  =  06. 

The  available  data  for  evaluating "the  Fc(k),  Gc(k),  Fc^(k), 

and  Gc^(k)  functions  have  been  given  in  various  forms,  beginning  with 

the  work  of  Possio.  One  form  in  which  all  of  the  available  data  can  be 
put  is  as  follows: 

Translatory  motion: 


L(s)  =  jrpcV2eiks 


|(zi  +  1Z2) 


M(s)  =  Jtpc2V2elkS  +  iM2) 


Pitching  motion  (including  contribution  due  to  pitching  velocity): 


L(s)  =  -JtpcV2elks  ©(Z3  +  iZ^) 


M(s)  =  -Jtpc2V2elks  e(^M3  +  iM^ 


where  the  moment  is  taken  about  the  quarter-chord  point  and  the  axis'  of 
rotation  for  the  harmonically  pitching  airfoil  is  located  at  the  leading 
edge . 

In  order  to  make  use  of  the  reciprocal  equations  (5a)  and  (5b)  to 
determine  the  indicia!  lift  function  due  to  pitching  velocity  kig(s), 

the  component  of  the  flutter  coefficients  due  to  pitching  velocity  must 
first  be  separated  from  the  flutter  coefficients  given  by  equations  (8a) 
and  (8b)  and  then  converted  to  the  form  given  by  equation  (6).  If  the 
expression  for  h  determined  by  equation  (4)  is  substituted  into 
equation  (7a)  and  if  the  resulting  expression  is  subtracted  from  equa¬ 
tion  (8a),  the  following  expression  for  the  lift  can  be  obtained  for  the 
component  due  to  pitching  velocity  on  a  harmonically  oscillating  airfoil: 


NACA  TN  2613 


Lq(s)  =  -JtpcV2elks  G  Z  -  ||  +  i(z,,  +  Zl 


(Z4  +  S)] 


4  t  oir;i  (9) 

Comparison  of  equations  (6)  and  (9)  leads  to  the  following  expressions 
for  the  functions  FCq(k)  and  GCq(k)  in  terms  of  the  oscillatory 

coefficients  Zlf  Z2,  Z3,  and  Z^ 


F~  (k) 


=  i(z4  +  2^) 


(10a) 


and 


%(k)  --s(Z3  -  ft) 


(10b) 


Substitution  of  equations  (10)  into  equations  (5)  then  allows  for  the 
determination  of  the  indie ial  function  kq  (s)  for  pitching  velocity. 

The  procedure  shown  for  determining  the  indicial  lift  function  due 
to  pitching  velocity  from  the  oscillatory  coefficients  may  be  applied 
to  determine  the  indicial  moment  function  due  to  pitching  velocity  mq  (s) 

defined  as  follows:  ^ 

M(s)  =  2jtpc2V2cp1  (s)  (11) 

where  q  is  the  magnitude  of  the  pitching  velocity.  The  reciprocal 
equations  can  be  written  for  the  moment  in  terms  of  the  oscillatory 
coefficients  and  the  indicial  moment  function  mi  (s)  as  follows: 


2  r00  Mq(k)  sin  ks 

~  k -  “ 


(s  >  0) 


(12a) 


and 


mi  (s)  =  Mq(0)  +  | 


Nq(k)  cos  ks 


0 


k 


dk 


(s  >0) 


(12b) 


where  Mq(k)  and  W^(k)  are,  respectively,  the  in-phase  and  out-of- 

phase  moment  coefficients  on  an  oscillating  airfoil  due  to  the  angular 
pitching  velocity  about  the  leading  edge  and  are  defined  by  the  equation 


10 


NACA  TN  2613 


Mn(s)  =  3tpc2V2eiks(2ik0)  |Mn(k)  +  iNn(k) 


In  equations  (12a)  and  ( 12b)  the  functions  M^(k)  and  N^(k)/k  must 
he  continuous  and  finite  in  the  interval  from  k  =  0  to  k  =  00  . 

In  a  manner  similar  to  that  shown  for  the  lift  case,  expressions 
for  Mq(k)  and  N^(k)  may  he  obtained  in  terms  of  the  oscillatory  * 

coefficients  M-^,  Mg*  Mg,  and  Mj^  as  follows: 


Mc/k)  "  “2k(M^  +  2k 


(14a) 


1  /  M2' 

Vk>  -  5k  M3  -  5k 


(llrt>) 


Summary  of  Available  Flutter  Coefficients  and  Related 
Data  for  a  Mach  Number  of  0.7 

The  two-dimensional  compressible  flutter,  coefficients  for  the  real 
and  imaginary  parts  of  the  lift  and  moment  for  sinking  and  pitching 
motion  are  given  in  table  I  for  a  Mach  number  of  0.7*  The  results  given 
in  table  I  were  taken  from  three  sources  (references  h  to  6)  and  have 
been  converted  to  the  form  given  by  equations  (7)  and  (8)  for  the  lift 
and  moment.  The  range  of  reduced  frequencies  taken  from  each  source  is 
indicated  in  table  I. 

Expressions  for  the  end  points  Fc(0)  and  Fc(°°)  can  be  deter¬ 
mined  independently  of  the  flutter  coefficients.  The  value  of  Fc(0) 
which  represents  the  steady-state  lift  is  determined  by  the  Prandtl- 

Glauert  factor  —  ^  (where  M  refers  to  the  Mach  number  of  the 

yrr w 

free  stream).  The  value  of  Fc(oo)  which  was  shown  in  reference  2  to 
correspond  to  the  value  of  k-j_(  s)  at  s  =  0  can  be  determined  by  the 
following  equation  given  in  reference  2  for  any  Mach  number: 
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The  corresponding  expression  for  the  moment  M(«>)  is  noted  in  refer- 
ence  2  as 


For  the  case  of  a  pitching  airfoil  the  expression  for  the  end  points 
for  the  Fc^(k)  and  Mq(k)  functions  may  he  obtained  in  a  manner  similar 

to  that  shown  for  the  case  of  sinking  airfoil.  For  an  airfoil  exper¬ 
iencing  a  sudden  change  in  pitching  velocity  about  its  leading  edge, 
the  steady-state  value  of  the  k]_^(s)  function  (which  can  be  shown  to 

correspond  to  FC(^(0))  can  be  determined  from  incompressible  flow 

together  with  a  correction  given  by  the  Prandtl-Glauert  factor.  From  the 
work  of  either  Theodorsen  or  Wagner  (reference  8  or  9)  the  incompres¬ 
sible  steady-state  value  can  be  determined,  and,  together  with  the 
Prandtl-Glauert  factor,  the  resulting  expression  for  FC(^(0)  for 

subsonic  flow  is  as  follows: 


(0)  =  knH  =  J 


1  -  M‘ 


Similarly,  for  the  moment  the  expression  for  Mq(0) 


Mq(0) 


,  ,  1  1 
V"'  ’•I5rr 


The  values  for  Fc^(oo)  and  Mq(oo)  may  be  obtained  from  the  starting 

values  of  the  k-^(s)  and  mq^( s)  functions,  respectively.  The  same 

method  of  analysis  for  obtaining  the  starting  values  of  the  k^s) 
and  m1(s)  functions  given  in  reference  9  for  a  sinking  airfoil  can  be 
applied  to  the  case  of  a  pitching  airfoil  to  obtain  the  starting  values 

for  the  k-|_  (s)  and  m^  (s)  functions  if  a  substitution  of  the  boundary 
Q.  Q. 

conditions  due  to  pitching  velocity  is  made  in  lieu  of  sinking  velocity 
as  indicated  by  equation  (l).  If  the  perturbation  velocity  component 
due  to  a  sudden  change  in  pitching  velocity  is  substituted  for  the 
perturbation  velocity  due  to  sinking  motion  in  the  analysis  given  in 
reference  9  for  determining  k^s)  and  mp(s)  at  s  =  0,  the  following 
expressions  are  obtained  for  the  starting  values  of  the  kp^s)  and 

mlq(s)  functions 


v°> 


RESULTS  AND  DISCUSSION 

Numerical  Solution  of  the  Reciprocal  Equations 
for  a  Mach  Number  of  0.7 


As  noted  in  reference  2,  the  indicial  lift  and  moment  functions 
determined  by  the  reciprocal  equations  containing  the  in-phase  lift  or 
moment  functions  (Fc(k)  or  M(k))  provided  a  more  reliable  and  simpler 
solution  than  the  functions  determined  by  the  reciprocal  equations  con¬ 
taining  the  out-of -phase  components  (Gc(k)  and  N(k)).  Consequently, 
the  indicial  lift  and  moment  functions  for  the  pitching-velocity  case 
were  determined  numerically  by  using  equations  (5a)  and  (12a)  in  a  man¬ 
ner  similar  to  that  shown  for  the  case  of  the  sinking  airfoil. 

In  figure  1  a  plot  is  shown  of  the  complex  oscillatory  lift  func¬ 
tions  due  to  pitching  velocity  at  a  Mach  number  of  0.7  which  were 
evaluated  from  the  flutter  coefficients  by  use  of  equations  (lOa)  and 
( 10b) .  If 

fq<k)  =  Fcq<k)  -  FCq(“) 


(23) 


NACA  TN  2613 


13 


is  substituted  into  equation  (5a),  the  following  form  which  can  be 
evaluated  graphically  more  readily  is  obtained: 


kiq(s)  =FCqH  .+  - 


°°fq(k)  sin  ks 


dk 


(24) 


A  plot  of  .the  function  fq(k)  is  shown  in  figure  2.  The  integrand  in 
equation  (24)  was  graphically  evaluated  for  several  values  of  the  parame 
ter  s.  In  figure  3>  plots  of  the  indicial  lift  function  k^  (s)  are 

shown  for  Mach  numbers  of  0.7  and  0.  The  indicial  lift  function  k^s) 
due  to  angular  position  alone  (corresponding  to  the  case  of  a  sinking 
airfoil)  is  plotted  in  figure  4  as  determined  from  reference  2  for  a 
Mach  number  of  0.7-  Also  shown  in  this  figure  are  the  beginning  posi¬ 
tion  of  the  curve  for  a  Mach  number  of  0.7  as  determined  by  the  equa¬ 
tions  given  in  reference  the  curve  for  a  Mach  number  of  0.8  given 
in  reference  3>  and  the  curve  for  a  Mach  number  of  0.  Comparison  of 
the  k]_(s)  and  k^(s)  functions  at  Mach  numbers  of  0  and  0.7  indi¬ 
cates  that  the  growth  of  lift  to  the  steady  state  for  a  pitching  airfoil 
is  less  rapid  in  subsonic  compressible  flow,  than  for  incompressible  flow 


The  indicial  moment  function  mi^(s)  due  to  a  sudden  change  in 

pitching  velocity  may  be  obtained  in  a  manner  similar  to  that  shown 
for  determining  the  indicial  lift  function  kiq(s).  In  figure  5,  plots 

are  shown  of  the  complex  moment  function  due  to  pitching  velocity  for  a 
Mach  number  of  0.7  determined  from  the  flutter  coefficients  M]_,  Mg, 
M3,  and  M^  by  equations  (l4).  The  asymptotic  value  for  Mq(k)  is 
determined  by  equation  (22) .  Because  of  the  nature  of  the  curves  at 
the  higher  reduced  frequencies,  the  curves  are  represented  by  dashed 
lines  to  indicate  a  degree  of  unreliability.  Reference  2,  however, 
indicates  that  errors  in  the  curves  at  .the  higher  reduced  frequencies 
would  have  little  effect  on  the  indicial  function. , 

In  a  manner  similar  to  that  shown  for  the  lift  case,  the  M«(k) 
function  was  transformed  by  the  following  equation 

mq(k)  =  Mq(k)  -  Mq(»)  (25) 

in'  order  to  facilitate  the  numerical  determination  of  the  indicial 
moment  function  mx  (s').  A  plot  of  mq(k)  is  shown  in  figure  6.  The  , 

function  miq(s)  was  evaluated  graphically  by  substitution  of  equa¬ 
tion  (25)  into  equation  (l2a)  for  several  values  of  the  parameters  and 
is  plotted  in  figure  7  together  with  the  solution  for  a  Mach  number 
of  0.  In  figure  8  a  plot  of  the  indicial  moment  function  due  to 


Ik 
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angle -of- attack  change  (no . pitching  velocity)  mp(s)  is  shown  for  Mach 
numbers  of  0.7  and  0.8  determined,  respectively,  from  references  2 
and  3  together  with  the  solution  for  a  Mach  number  of  0.  - 

References  8  and  9  show  that  the  indicial  lift  due  to  pitching 
velocity  on  an  airfoil  pitching  about  its  three -quarter- chord  position 
is  impulsive  at  s  =  0  and  is  zero  thereafter  for  incompressible  flow. 
For  compressible  flow  the  indicial  lift,  however,  is  initially  finite 
and  is  time -dependent  thereafter,  as  is  shown  subsequently.  When  the 
results  presented. in  figures  3;  ^  1,  and  8  are  used,  the  indicial  lift 
and  moment  functions  due  to  pitching  velocity  for  an  airfoil  pitching 
about  any  axis  can  be  obtained  readily  by  the  following  general  trans¬ 
formations: 

(klq)  (s)  =  klg( s)  -  “  kl(s)  (26a) 

and 

(mih  (s)  =  ni!  (s)  -  ^  m^s)  (26b) 

\  4/x  q  c' 

where  x  denotes  the  location  of  the  axis  of  rotation  measured  posi¬ 
tively  to  the  right  from  the  leading  edge  of  the  airfoil  (see  sketch  in 

section  entitled  "Method  of  Analysis") .  For  the  case  of  an  airfoil 

pitching  about  its  three -quarter- chord  location,  equations  (26a)  and 
(26b)  become  , 

(klq)3£(s)  =  klq(s)  kl(s)  (27a) 

'  Ik 

and 

•/mi  \  (s)  =  mlq(s)  -  m^s)  (27b) 

't 


The  functions  (ki  \  (s)  and  /mi  \  (s)  were  determined  for  a 

\  V3c  \  <i)3c 

Mach  number  of  0*7  from  equations  (27a)  and  (27b)  and  from  the  results 
plotted  in  figures  3>  1>  and  8.  Plots  of  the  functions  fk]\  (s) 

1  q% 
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and  (mi  )  (s)  are  shown  in  figures  9  and  10,  respectively,  together 

V  q/3£ 

4 

with  the  'indicial  functions  for  the  case  for  a  Mach  number  of  0. 

Also  included  in  these  figures  are  the  solutions  for  a  Mach  number 
of  0.8  and  part  of  the  solution  for  a  Mach  number  of  0.7,  both  deter¬ 
mined  from  reference  3-  A  comparison  of  the  (s)  function  for 

4 

subsonic  compressible  and  incompressible  flow  indicates  that,  although 
this  function  is  impulsive  at  s  =  0  and  zero  for  s  >  0  for  incom¬ 
pressible  flow,  the  function  is  finite  at  s  =  0  and  decays  to  zero 
very  rapidly  in  the  compressible  case.  If  for  compressible  flow  the 
circulation  is  assumed  to  be  zero  for  s  >0  (as  is  the  case  for  this 
component  at  M  =  0) ,  the  time-dependent  function  which  is  present  in 
compressible  flow  may  be  attributed  to  the  time -dependent  apparent-mass 
effects.  Also,  the  part  of  the  k^  (s)  function  shown  in  figure  4 

for  s  greater  than  approximately  4  may  be  associated  only  with  the 
lift  due  to  circulation.  Comparison  of  this  part  of  the  curve  with  the 
curve  for  a  Mach  number  of  0  in  the  region  s  >  4  indicates  that  the 
growth  of  lift  to  the  steady  state  is  relatively  less  rapid  for  com¬ 
pressible  flow  than  for  incompressible  flow. 


Approximation  of  Indicial  Lift  and  Moment  Functions 
by  Analytic  Expressions 

Since  the  exponential  function  has  a  simple  operational  equivalent 
and  it  has  been  found  convenient  to  approximate  the  ki(s)  function  at 
a  Mach  number  of  0  (see  reference  10) ,  a  limited  series  of  such  func¬ 
tions  were  chosen  to  approximate  the  indicial  lift  and  moment  functions 
for  sinking  and  pitching  motion  at  a  Mach  number  of  0.7.  The  functions 
were  found  to  fit  these  curves  quite  well  and  are: 

kj_(s)-  =  1.4(l  -  0.364e'0,0536s  -  0.405e'0,357s  +  0.4l9e"°-9°2s)  (28) 


k-,  \  (s)  =  -0.083e 

<l/3c 


-°-800s  -  O^Se"1'*56*8  +  0.l49e-2-^s 


(s)  =  -0.2425e-°-974s  +  0.084e-°-668s  -  O^e-0-4385 


(30) 
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/mlq)3c(s)  -  -0.0875(l  ♦  O.llltle-0-181555  -  1.233e'1-1,'1S  *  O.SSSTe'1*-0^) 

T 

(3i) 

These  analytic  approximations  given  by  equations  (28),  (29),  (30), 
and  (31)  are  plotted  in  figures  k,  9,  8,  and  10,  respectively.  Com- 
parison  of  the  approximation  with  the  actual  curves  in  these  figures 
indicates  the  order  of  good  agreement  reached. 

The  corresponding  approximate  expressions  for  the  harmonically 
oscillating  airfoil  can  he  found  from  the  reciprocal  form  of  equa-  - 
tions  (2)  and  (5)  where  the  harmonically  oscillating  functions  are 
expressed  in  terms  of  their : corresponding  indicial  function.  For  the 
case  of  the  sinking  airfoil,  the  reciprocal  form  of  equation  (2)  is 

P°°  -iks 

Fc(k)  ■+  iGc(k)  =  ik  /  -k^sje  ds  (32) 

J0 


If  ik  is  considered  as  the  operator  in  the  Laplace  transformation, 
then  equation  (32)  is  simply  ik  times  the  Laplace  transformation 
of  kp(s).  Therefore  Fc(k)  and  Gc(k)  are,  respectively,  the  real 
and  imaginary  parts  of  ik  times  the  Laplace  transform  of  k^(s). 
Substitution  of  equation  (28)  into  equation  (32)  leads  to  the  following 
expressions  for  Fc(k)  and  Gc(k)  : 


Fc(k)  =  1.4 


1  - 


0 . 364k5 


0 . 405k2 


0.4l9k2 


(0.0536)2  +  k2  (0.357) 2  +  k2  (0.902)2  +  k‘ 


>  (33) 


Gc(k)  =  1.4 


(0- 364)  (0. 0536) -k 

(0.0536) 2  +  k2 


(0.405) (0. 357) k  ,  (0.419) (0.9Q2)k 

7oT357)2~+^  (0.902)2  +  k2JJ 


In  a  manner  similar  to  that  shown  for  the  lift  case,  the  following 
approximate  expressions  for  M(k)  and  N(k)  can  he  obtained: 
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M(t)  =  .  +  o.oBta2 


0.069k2 


(0.974) 2  +  k2  (  0.668)2  +  k2  (0.438)2  +  k2 


,,  v  _  ( -0.2429) (0. 974) k  (0.084) (0.668)k  _  (0.069) (o.438)k 
(0.974) 2  +  k2  (0.668) 2  +  k2  (0.438)2  +  k2 


(34) 


j 


The  expressions  for  the  /f_  ]  (k)  and  (Gr  \ 

V  <jj3c  lC<3/3c 

4 


(k)  lift  functions 


due  to  pitching  velocity  are 


N 


(k)  = 


0.083k‘ 


0.293k2  + _ o.i49k2 


3c  (0.800)2  +  k2  (I.565)2  +  k2  (2.44)2  +  k2 
4 


N 


(k)  = 


3c 

4 


(-0.083) (0.800)k  _  (0.293) (l-969)k  ^  (0.l49)(2.44)k 
(0.800)2  +  k2 


(I.565)2  +  k2  (2 . 44) 2  +  k2 


>  (35) 


and  the  expressions  for  the  ^Mqj^c(k)  and  ^Nqj^c(k)  moment  functions 


are 


^MqJ3c(k)  =  -0.0875 1 1  + 

4 


0.1l4lk‘ 


1.233k2  +  013337k2 


3c  I  (0.1865)2  +  k2  ( 1 . l4l) 2  +  k2  (4.04)2  +  k2 


(”4<k)  ■  '°-0875fe 

4 


Il4l)(0.l865)k  (l.233)(l.l4l)k 


1865) 2  +  k2  (1.141)*  +  k; 


2  .  .,.2 


(0-3337)  (4. 04)k 
(4.o4)2  +  k2 
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CONCLUDING  REMARKS 

The  indicial  lift  and  moment  functions  due  to  a  sudden  change  in 
pitching  velocity  have  been  obtained  from  the  available  data  on  a 
harmonically  pitching  airfoil  at  a  Mach  number  of  0.7  by  the  use  of 
reciprocal  relations.  Comparison  of  the  results  obtained  for  the  case 
of  a  Mach  number  of  0-7  with  incompressible  flow  indicates  that  the 
growth  of  lift  to  the  steady  state  appears  to  be  less  rapid  for  com¬ 
pressible  flow  than  for  incompressible  flow  for  a  pitching  airfoil. 
Although  the  circulatory  lift  component  due  to  pitching  velocity  for 
an  airfoil  rotating  about  its  three -quarter- chord  point  was  zero  for 
incompressible  flow,  a  time -dependent  component  of  lift  was  found  to 
exist  for  compressible  flow  and  decayed  to  a  negligible  value  in 
approximately  2  chords. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  October  30,  1953 
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Figure  4.-  Comparison  of  indicial  lift  functions  due  to  a  sudden  change 
in  angle  of  attack  (without  pitching  motion)  at  M  =  0,  M  =  r 
M  =  0.8. 
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